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ABSTRACT 

We propose that the anomalously bright white dwarf luminosity function observed in NGC 6791 
(Bedin et al 2005) is the consequence of the formation of 0.5M Q white dwarfs with Helium cores instead 
of Carbon cores. This may happen if mass loss during the ascent of the Red Giant Branch is strong 
enough to prevent a star from reaching the Helium flash. Such a model can explain the slower white 
dwarf cooling (relative to standard models) and fits naturally with scenarios advanced to explain Extreme 
Horizontal Branch stars, a population of which are also found in this cluster. 

Subject headings: open clusters and associations: individual (NGC 6791) - stars: white dwarfs, 
luminosity function, mass loss, horizontal branch, evolution 



1. INTRODUCTION 

The open cluster NGC 6791 is an interesting object for 
several reasons. It is amongst the richest of the Galactic 
disk clusters, possibly the oldest, and possesses a signifi- 
cantly supersolar metallicity. These attributes have made 
it a target of several detailed studies. Nevertheless, there 
is still some uncertainty as to the true age. Stetson et al 
(2004) give an estimate ~ 12 Gyr, while studies such as 
that of Carney, Lee & Dodson (2005) or King et al (2005) 
yield estimates closer to ~ 8 Gyr. The age uncertainties 
are covariant with uncertainties in the cluster metallicity 
([Fe/H] - 0.25-0.5) and distance. 

The above estimates are all arrived at by measure- 
ments of the main sequence turnoff and giant sequences. 
The detection of a significant white dwarf population in 
NGC 6791 (Bedin et al 2005) is thus of considerable in- 
terest, since it offers the potential to measure cluster pa- 
rameters by entirely independent means. However, fitting 
standard white dwarf models to the observed luminosity 
function yields a cluster age of ~ 3 Gyr. In this paper 
we examine some non-standard white dwarf models and 
propose an explanation for this puzzling discrepancy. 

In § [21 we review the results of Bedin et al and sum- 
marise the various possible explanations they were able to 
rule out. Thereafter we describe two additional possibili- 
ties not considered in that paper (residual nuclear burning 
and the production of massive white dwarfs with Helium 
cores) . In § we compare these models to the data and 
demonstrate that the latter may indeed provide an expla- 
nation for the NGC 6791 luminosity function. In §0]we 
examine some of the consequences of this model and pos- 
sible predictions. 

2. THE WHITE DWARF LUMINOSITY FUNCTION 

The luminosity function derived by Bedin et al has the 
characteristic shape expected for the white dwarf popu- 
lation drawn from a burst of star formation - strongly 
peaked near a limiting value followed by a sharp drop at 
lower luminosities. The problem is that the location of 
the peak is far brighter than was anticipated. The ex- 
pected absolute magnitude for a O.5M white dwarf of 
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age 8 Gyr is Mqo§ ~ 16 which, combined with a distance 
modulus of /i606 — 13.44, yields an expected peak loca- 
tion of -F6061U = 29.5. Yet the luminosity function of 
Bedin et al peaks at approximately F606W ~ 27.5. They 
present compelling evidence that the peak lies well above 
their completeness limit, so that one is forced to consider 
ways to generate a peak in the luminosity function at much 
brighter magnitudes than expected. 

Bedin et al reviewed several possible explanations for 
their luminosity function. Simply decreasing the distance 
or extinction is not viable, as it would drive the age de- 
rived from the main sequence turnoff to unacceptably large 
values. Small changes in the white dwarf models, such as 
changing the Hydrogen layer thickness or the ratio of Car- 
bon and Oxygen in the core, do not produce a big enough 
effect. White dwarfs produced by truncated stellar evo- 
lution in binaries (Kippenhahn, Kohl & Weigert 1967) do 
cool more slowly, but generally have a range of masses 
<~ 0.3-0.4M Q and so would be redder than the observed 
white dwarf cooling track. Indeed, any explanation in- 
volving binary stars needs to account for the narrow mass 
range and the lack of bright (relative to a white dwarf) 
companions at the present time. 

There is one potential explanation in the literature, that 
actually predates the Bedin et al observation. Bildsten & 
Hall (2001) and Deloye & Bildsten (2002) discuss the re- 
tardation of white dwarf cooling that is made possible by 
the sedimentation of 22 Ne during the cooling process. Al- 
though the contribution is small for most white dwarfs, 
it gets larger if there is more 22 Ne - which is expected 
from more metal-rich systems. Thus, Deloye & Bildsten 
in fact predicted that the effects of sedimentation would be 
largest in such a system as NGC 6791. At first glance, the 
NGC 6791 white dwarfs might appear to be a stunning 
confirmation of the prediction of 'sedimentars', but the 
size of the predicted effect is somewhat smaller than that 
needed to explain the Bedin et al result. Figure 12 and 
13 of Deloye & Bildsten can be converted into a prediction 
that the bulk of the white dwarfs should be found between 
F606W = 28-29, depending on the value of the assumed 
diffusion coefficient (the brightest value being found for 



2 



NGC6791 White Dwarfs 



a rate ten times faster than the nominal but uncertain 
value). Nevertheless, the agreement may be improved by 
further calculations in progress (L.Bildsten, personal com- 
munication) . 

At present, however, the observations still require an ex- 
planation and so we wish to re-examine two issues touched 
on by Bedin et al, but not explored to the fullest. 

2.1. Thick Hydrogen Layers and Nuclear Burning 

Bcdin et al did examine the effects of changing the Hy- 
drogen layer mass on the cooling, but as far as can be 
told from the paper, considered only the effect it has on 
heat transport from the core to the surface. For a large 
enough Hydrogen layer, the pressure and temperature at 
the base are high enough to maintain some level of resid- 
ual nuclear burning by the pp process, i.e. an additional 
heat source. This mechanism has been invoked in previ- 
ous instances of anomalously young white dwarfs with low 
mass (e.g. Alberts et al 1996; Driebe et al 1998). Could 
residual nuclear burning be the explanation for the bright 
white dwarfs in NGC 6791? Unfortunately, this seems 
unlikely. The effects of nuclear burning are only signifi- 
cant for masses M < 0.25A/ Q (Driebe et al 1999; Althaus, 
Serenelli & Benvenuto 2001; Hansen, Kalogera & Rasio 
2003). 

To confirm this, we calculate several new white dwarf 
cooling sequences using the code described in Hansen 
(1999). The Helium layer is taken to be qne = 10~ 2 . 
All models have total mass 0.5M Q . We use the same C/O 
profile as in the previous models, but consider a variety 
of Hydrogen layer thicknesses, ranging from qn = 10~ 4 to 
qn = 5 x 10~ 3 . The lower end of the range is the canonical 
value expected for most white dwarfs, but larger values are 
possible given the uncertain nature of the mass loss history 
of evolved stars. 

In the initial stages (central temperatures > 10 s K) 
of the approach to the cooling sequence, the evolution is 
driven by neutrino cooling in the core and the cooling is 
similar for all models. However, as the star shrinks towards 
a more compact, degenerate configuration, the pressure 
and temperature at the base of the Hydrogen layer increase 
and nuclear burning is possible for some range of Hydrogen 
layer masses. For even the most massive layers considered, 
the models remain brighter than the fiducial (qjj — 10 -4 ) 
models only for ~ 2 Gyr. This is because, if one increases 
the mass of the Hydrogen layer, the pressure and density 
at the base of the envelope are higher and so the rate of 
burning is higher. In the end the effect of a larger reser- 
voir is balanced by a higher rate of consumption, limiting 
the delay one can achieve. After 2 Gyr, no model has 
a Hydrogen layer larger than qn ~ 5 x 10~ 4 . Thus, af- 
ter ~ 8 Gyr, the white dwarf will have approximately the 
same luminosity regardless of the initial Hydrogen mass. 
This is demonstrated in Figure^ The white dwarf cooling 
time (neglecting the main sequence lifetime of the progen- 
itor for now) for white dwarfs at .F606IU ~ 27.3 is 1.6 Gyr 
for a standard C/O model. Even the thickest Hydrogen 
envelope considered (qjj ~ 5 x 10~ 3 ) only lengthens the 
cooling time to this luminosity by 0.6 Gyr. 

2.2. Helium Cores 

The notion that the NGC 6791 white dwarfs have cores 
composed of Helium, rather than Carbon or Oxygen, was 



touched on briefly by Bedin et al. However, they ruled 
this out as Helium core white dwarfs are believed to result 
from binary evolution and the majority have masses that 
range ~ 0.3-0.4M Q and would thus be too red to fit the 
observations. 

We believe this to be an overly conservative restriction. 
In this section we examine models for white dwarfs of mass 
O.45-O.55M0 (models with this mass are consistent with 
the colour and magnitude of the upper cooling sequence) 
which nevertheless possess cores composed purely of He- 
lium. We will leave issues of provenance to §0 

For white dwarfs hot enough to be far from the strongly 
coupled regime in the core, the heat content (and so the 
cooling time) of the star is inversely proportional to the 
mass number of the core constituent. Thus, if we consider 
two 0.5M Q white dwarfs, one composed of Carbon and the 
other of Helium, both with an age 8 Gyr, then the latter 
will be considerably brighter. In fact, it will have approxi- 
mately the same luminosity as the Carbon core model has 
at an age ~ 8/3 ~ 2.7 Gyr. This is approximately the 
age inferred by Bedin et al. Figure Q shows that the He- 
lium core model fares considerably better than the C/O 
models, yielding a cooling time of ~ 4 Gyr to My ~ 27.3. 
While this is still somewhat below the 8 Gyr age of the 
cluster, we must recall that the cluster age is the sum of 
the time spent on both the cooling sequence and on the 
main sequence. Thus, the difference between the cluster 
age and the Helium-core cooling time indicates the mass 
of the main sequence progenitors of this anomalous popu- 
lation. 

We use the analytic stellar evolution formulae encoded 
in the SSE package (Hurley, Pols & Tout 2000) to calculate 
the main sequence lifetimes of stars with Z = 0.035. For 
the low end of the cluster age range (~ 8 Gyr), the pro- 
genitor mass for Helium core white dwarfs at the observed 
cutoff is ~ 1.45M Q , while for the high end of the age range 
(~ 12 Gyr), the progenitor mass is ~ 1.17M Q . Exploring 
the range of applicable white dwarf masses (O.45M to 
0.55M Q ) does not change the conclusions significantly. 

Next we will try to model the luminosity function. Be- 
fore doing this we need to consider in more detail the pos- 
sible origins of such a population. 

3. AVOIDING THE HELIUM FLASH IN NGC 6791 

There are good reasons to believe most white dwarfs 
have C/O cores. Upon reaching the tip of the Red Gi- 
ant Branch (RGB), a star ignites core Helium burning 
(under degenerate conditions for lower mass progenitors - 
which leads to a thermonuclear runaway, a.k.a. the Helium 
flash), and the star moves onto the Horizontal Branch. It 
undergoes extended core Helium burning followed by shell 
burning on the Asymptotic Giant Branch before becom- 
ing a white dwarf. In order to avoid converting the Helium 
core to Carbon and Oxygen the star has to find some way 
to circumvent these later stages of stellar evolution. 

Stars in binaries sometimes manage this feat by losing 
their envelope on the ascent of the RGB due to Roche lobe 
overflow i.e., mass transfer to a companion. However, this 
leads to a range of masses, most of which are consider- 
ably lower than the - O.5M needed to fit the NGC 6791 
cooling sequence. Bedin et al have already considered and 
rejected this possibility. In order for our model to work we 
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have to hypothesize that many single stars in NGC 6791 
have managed to lose their envelopes before reaching the 
Helium flash, so that they move directly from the RGB to 
the white dwarf cooling sequence. 

In fact, this hypothesis dovetails quite naturally with 
several theories for the origins and nature of Extreme Hor- 
izontal Branch (EHB) stars. This term is used to de- 
scribe a class of stars, found both in the field and in clus- 
ters, which appear to be related to traditional Horizon- 
tal Branch stars but which are hotter/bluer. It has been 
suggested (Faulkner 1972; Sweigart, Mengel & Demarque 
1974) that such stars are core Helium burning stars with 
particularly thin Hydrogen envelopes, possibly as a result 
of extreme mass loss on the RGB. Subsequent work has 
developed this picture even further. Castellani & Castel- 
lani (1993) and Castellani, Luridiana & Romanicllo (1994) 
report the formation of what they call "Red Giant Strag- 
glers" in models for globular cluster evolution. In some 
cases, the stars make it to the white dwarf cooling se- 
quence before igniting Helium. D'Cruz et al (1996) con- 
sidered models with a range of mass loss rates on the RGB 
and found that, in addition to the formation of EHB stars, 
they formed so-called "flash-manque" stars, which lose so 
much mass that they never ignite Helium and simply go di- 
rectly to Helium core white dwarfs. In fact, given the large 
amount of mass loss required to form EHB stars, especially 
the very bluest 'blue hook stars' (Brown et al 2001; Cassisi 
et al 2003; Moehler et al 2004), it would require extreme 
fine tuning to produce EHB stars without also producing 
sonic Helium core white dwarfs. This is important because 
NGC 6791 possesses a significant population of such EHB 
stars (Kaluzny & Udalski 1992; Liebert, Saffer & Green 
1994). Thus, if one takes the above models at face value, 
a substantial population of Helium core white dwarfs is 
expected wherever EHB stars are found and thus should 
be found in NGC 6791. 

We now try to model the observed luminosity function 
using this framework. We note that not all the white 
dwarfs in NGC 6791 can have Helium cores. The models 
discussed above that successfully avoid the Helium flash 
do so with progenitors that begin their lives with masses 
slightly larger than 1M . Thus, in our models we shall 
impose a critical mass m cr it, above which stars always 
produce standard C/O core white dwarfs. Furthermore, 
NGC 6791 does possess normal Helium-burning stars (the 
EHB stars make up ~ 15% of the Helium-burning stars ac- 
cording to Liebert et al 1994) and even the EHB stars are 
Helium-burning, so that clearly some C/O white dwarfs 
are being produced. It seems likely that stellar evolution 
in this cluster explores all three post-RGB avenues dis- 
cussed above. One may estimate the branching ratios as 
follows. 

King et al (2005) note that two of the EHB candi- 
dates from Kaluzny & Udalski (1992) lie within the HST 
field. Models suggest that this evolutionary stage lasts for 
~ 10 8 years. The offspring of EHB stars are C/O white 
dwarfs and so, using the cooling time of C/O models (to 
F606W = 28), we estimate that - 2 x 10 9 /10 8 x 2 = 40 of 
the white dwarfs in this field brighter than F606W = 28 
should have come through this channel. In addition, C/O 
white dwarfs are also the end product of normal Helium- 
burning stars in the cluster. Although we do not have a 
strict count of HB stars in this field, Liebert et al (1994) es- 



timate that ~ 15% of all NGC 6791 Helium-burning stars 
are EHB stars. Thus, we estimate that the total number of 
C/O white dwarfs produced in this field with F606W < 28 
is ~ 40/0.15 ~ 270. The total number of white dwarfs ob- 
served by Bedin et al in the same field is ~ 600. The 
difference between these two numbers is the number of 
Helium core dwarfs that went directly from RGB to cool- 
ing sequence. Thus, the ratio of C/O core white dwarfs to 
Helium core dwarfs is estimated to be ~ 270 : 330 ~ 5 : 6. 
If we break this up further into EHB/normal HB/direct 
He core we estimate ratios ~ 1 : 6 : 8. The principal un- 
certainty in this procedure is the degree to which internal 
cluster dynamical evolution violates the implicit closed box 
assumption. Although EHB stars and white dwarfs should 
have very similar masses, the normal HB stars may have 
slightly higher masses and thus may be slightly underep- 
resented (since the field is away from the cluster core) in 
the census just described. In light of this uncertainty, we 
adopt a ratio of C/O to Helium cores of 1:1 below. 

Figure [3 shows the luminosity function of Bedin et al 
(now binned in 0.5 magnitude bins) compared to sample 
8 Gyr model luminosity functions. The best fit value of 
merit = 1.6M is slightly larger than the 1.45-M© quoted 
in § l2.2l because the peak of the realised Monte Carlo num- 
ber distribution lies slightly below the cutoff. The solid 
line shows a histogram for a model in which all stars with 
m < m cr it form Helium core white dwarfs. The dotted line 
shows the C/O core luminosity function for the same pa- 
rameters. The dashed line corresponds to a model in which 
50% of the stars with m < m cr i t form Helium cores and 
the rest form C/O cores. We have included the effects of 
incompleteness, using the results from Bedin et al, kindly 
provided by Ivan King. The truncation of the C/O core 
luminosity function shows where this becomes important. 
The truncation of the Helium core luminosity function oc- 
curs at significantly brighter magnitudes. The uncertain 
age of the cluster means that there is some flexibility in 
the choice of parameters. Figure [3] shows a fit using an 
age of 12 Gyr and m cr i t — 1.25M Q . The solid, dotted 
and dashed lines once again indicate luminosity functions 
for Helium core, C/O core and a 1:1 combination, as in 
Figure |3 

There are several points to note about these fits. We 
have not performed a detailed parameter scan in fitting 
the data because we do not possess a proper photometric 
uncertainty map for this data (which indicates the proba- 
bilistic relation between intrinsic and observed magnitude 
as a function of model magnitude - see Hansen et al (2004) 
for a more detailed discussion of this question in the con- 
text of the globular cluster M4). The structure in the 
CMD in Figure 1 of Bedin et al suggests this is likely to 
be an important consideration in performing a true fit. 
We also need to assume something about the distribu- 
tion of progenitor masses. We have assumed top heavy 
mass functions, with a slope x < — 1 (where Salpeter is 
x — 1.35). This is necessary to obtain luminosity func- 
tions as peaked as those seen here, and is also broadly 
consistent with the main sequence mass function observed 
by King et al (2005) for this cluster (the number of main 
sequence stars per magnitude bin increases with decreas- 
ing magnitude/increasing mass near the turnoff). For an 
open cluster as old as NGC 6791, there has likely been sig- 
nificant dynamical evolution and so one should be cautious 
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when interpreting this as a true mass function. 

4. DISCUSSION 

In this paper we propose a model to explain the origin of 
the peak in the NGC 6791 luminosity function at luminosi- 
ties somewhat brighter than those expected on the basis 
of traditional white dwarf cooling models. Our proposal is 
that these are the result of the formation of Helium core 
white dwarfs resulting from strong mass loss on the RGB, 
thereby avoiding an episode of core Helium burning on 
the Horizontal Branch. In support of this model we point 
to the substantial population of EHB stars in NGC 6791, 
which are proposed to be stars which lost enough mass 
to almost make it to the white dwarf cooling track with 
their Helium cores intact but finally ignited Helium burn- 
ing during the final contraction stage. In many models, 
EHB stars and Helium core white dwarfs are expected to 
form together as they are consequences of the same phe- 
nomenon - strong mass loss. 

The scenario discussed here does make a definite pre- 
diction - that the white dwarf luminosity function for 
NGC 6791 should be bimodal, with a second peak at 
fainter magnitudes resulting from traditional C/O white 
dwarfs which did, indeed, pass through the Horizontal 
Branch phase. Figure 0] shows the expected luminosity 
function for our simple 8 Gyr model with m cr it = 1.6Mq. 
The solid, dotted and dashed histograms show the same 
luminosity functions as in Figure [5] but now we have not 
modelled the effects of incompleteness. We see that the 
second peak is expected to lie at magnitudes F606 ~ 29- 
30. The exact value will vary depending on the mass of 
the white dwarfs at the faint end. The value shown here 
is probably a little optimistic since we used only 0.5M Q 
white dwarfs. 

The narrowness of the white dwarf cooling sequence sug- 
gests that, however these Helium-core white dwarfs form, 
it must be by a process that strongly favors core masses 
> OAMq. If heavy mass loss on the RGB is indeed the 
cause, this mass limit suggests that the progenitors must 



get within ~ 1 magnitude of the tip of the RGB before 
their evolution is truncated (based on metal-rich models 
using the formulae of Hurley ct al 2000). There is some 
observational support for this notion. Origlia et al (2002) 
used ISOCAM to search for infrared excesses suggestive of 
mass loss in the cores of five globular clusters. They found 
evidence for significant mass loss only close to the RGB 
tip. Similarly, Ita et al (2002) found many variable stars 
near the tip of the RGB in the Large Magellanic Cloud, 
possibly an indication of pulsation driven mass loss. While 
these results are encouraging, they only indicate that most 
of the mass loss on the RGB is indeed likely to occur near 
the tip. There is clearly still much work required to verify 
that the amount of mass loss on the NGC 6791 RGB is 
sufficient to justify our hypothesis. 

Finally, the referee has raised the question whether the 
existence of this anomalous white dwarf population poses a 
problem for the determination of population ages from the 
white dwarf cooling sequence (Fontaine, Brassard & Berg- 
eron 2001, Hansen 2004 and references therein). Within 
the model proposed here, this is not the case, because the 
difference between the Helium core and Carbon core white 
dwarfs is discrete in nature and any anomaly is readily ap- 
parent. Just as someone cannot be "a little bit pregnant", 
a star must belong to one population or the other. If 
the model of Deloye & Bildsten is correct, however, then 
the anomaly is "tunable" via the abundance of 22 Ne and 
thereby does indeed introduce an additional uncertainty 
to take into account. Thus, it is of considerable interest to 
distinguish between these two possibilities with a deeper 
observation of this cluster. 
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Fig. 1. — The horizontal dotted lines bracket the luminosity range corresponding to the Bedin et al peak. We have used the distance 
modulus Meo6 = 13.44 and colour— T e g transformations from Bergeron, Wesemael & Beauchamp (1995). The dot-dashed curve corresponds 
to a standard 0.5Mq white dwarf with C/O core and normal Helium and Hydrogen layers {que = 10 — 2 , qu = 10 — 4 ). The dashed curve is 
for the same model but now with a much thicker Hydrogen layer qu = 5 X 10 — 3 . This results in a period of nuclear burning during which 
the star is considerably brighter. However, the Hydrogen fuel is exhausted after ~2 Gyr. The solid curve shows a 0.5Mq model with pure 
Helium core and a Hydrogen layer envelope of qn = 10 — 4 . The cooling is ~ 2.5 times slower because of the higher total heat capacity of the 
Helium core. 

Fig. 2. — The points show the observed luminosity function from Bedin et al. (2005). This is compared to three theoretical luminosity 
functions. All are calculated for a cluster age of 8 Gyr. The solid luminosity function represents a population for which all stars with 
m < m cr it = I.GMq for Helium core white dwarfs. The dotted luminosity function is the standard model in which all stars produce C/O core 
white dwarfs. The dashed luminosity function is one in which only 50% of the stars with m < m cr it make Helium core white dwarfs. The 
theoretical models include a correction for the incompleteness of the Bedin et al observations. For simplicity, in this model all white dwarfs 
are assumed to have mass 0.5Mq. 

Fig. 3. — The format for this model is the same as for Figure|5] except that the cluster age is now assumed to be 12 Gyr and the value of 
m crit = 1.25M . 

Fig. 4. — The format for this figure is once again that of Figure|5] including the age and m cr it, except that now we make no correction for 
incompleteness. This shows the large number of white dwarfs we expect to find at fainter magnitudes if this model is correct. 
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